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Abstract 

The distribution function of the end-to-end distance of a semifiexible polymer, 
G{R\ L) (where R denotes the end-to-end distance and L the contour length), 
is calculated using a meanfield-like approach. The theory yields an extremely 
simple expression for G{R; L) which is in excellent agreement with Monte 
Carlo simulations. The second and fourth moments of G{R; L) agree with 
exact results for a semifiexible polymer in both the random coil and the rod 
limit. 
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Many properties of isolated, flexible polymer molecules are now understood |T[. For 
example, theoretical methods can be used to calculate the distribution function of flexible 
polymer molecules to very high accuracy. Unfortunately, many polymer molecules have too 
much internal stiffness to be successfully modeled as flexible chains 0. This is especially true 
for several important biopolymers such as actin, DNA, and microtubules [Q. A measure of 
the stiffness of a polymer is the persistence length, Ip. To understand the persistence length, 
consider two points on the polymer separated by a length / along the contour of the polymer 
backbone and construct the tangent vectors to this contour. The persistence length is the 
length along the contour over which these tangent vectors become uncorrelated. Thus, for 
/ < Ip the tangent vectors would have significant correlation and for / > Ip the tangent 
vectors would have little correlation. Typical values of Ip for biopolymers range from several 
nm to a few mm. If the contour length L of a polymer is of the same order of magnitude as 
Ip (or smaller) then the flexible chain model is inadequate for describing the polymer. For 
such a polymer it is imperative to include bending rigidity to describe the conformations of 
the chain. An appropriate model, called the semiflexible or wormlike model, was introduced 
in 1949 by Kratky and Porod [Q. The semiflexible model has been shown to provide a good 
starting point in the description of polymers with significant internal stiffness. 

In contrast to the situation of isolated, flexible polymers, there are many properties of 
semiflexible polymers that are not understood. Inspired by recent experiments that have 
probed the properties of semiflexible biological molecules, there has been a renewed interest 
in understanding their shapes. A central quantity for understanding the shape of a polymer 
is the distribution of end-to-end distances. Such a distribution function can, in principle, be 
measured by scattering experiments. Furthermore, it can be used to calculate the structure 
factor, which is needed as input in describing the dynamics of polymers. Therefore, finding 
a simple expression for this distribution function which is also accurate over a wide range of 
stiffness is important for further theoretical progress in understanding the properties of stiff 
chains. This paper provides just such a simple expression for the distribution of end-to-end 
distance of an ideal, semiflexible polymer chain. 

There have been several studies of the distribution of the end-to-end distance of semi- 
flexible chains Most of these studies have started by considering the chains near the 
rod limit and have computed corrections in powers of t^^, where t represents the ratio of the 
contour length to the persistence length. These calculations are very complicated and, more 
importantly, they do not provide reliable results in the interesting case where t is of the order 
one. Recently [^], Wilhelm and Frey [WF] have reported careful analytic and numerical 
(Monte Carlo) calculations for the radial distribution function of the end-to-end distance for 
a range of values of t. Their analytic expressions consist of infinite series of parabolic cylin- 
der functions (for two dimensions) and Hermite polynomials (for three dimensions). These 
expressions compare extremely well with their simulation results. Our approach differs from 
these previous studies. In this paper we use the method of our recent meanfield theory of 
semiflexible polymers to obtain a simple expression for the distribution function of the 



end-to-end distance. This expression is surprisingly accurate when compared to the analytic 
theory and simulations results for the range of t studied in [[l^ . 



The Hamiltonian for the semiflexible chain is taken to be [|T2| 
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where Ip (= n/ksT with k being the bending rigidity) is the persistence length, u(s) (= 
dr{s)/ds) is the unit tangent vector to the curve r(s) which describes the chain contour, s 
is the position measured along the chain contour, and L is the contour length of the entire 
chain. The distribution function of the end-to-end vector R is 



G(R; L) = m R - £ u{s)ds 



where the average is evaluated with respect to a thermal weight function, \l'[u(s)], 

^ JV[uis)]...nuis)] 
^■■■^ JV[u{smuis)] 

The thermal weight function for the semifiexible chain is 

^[u(s)] = 6{u{sf - 1) exp{~n[u{s)]) 
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In an earlier paper |jTT[ we showed that when the weight \E'[u(s)] is replaced by 



vI/mf[u(s)] = exp (^^j ds-\ Jjuisf - l]ds - e[u(0)^ - 1] - eHlf - 1]| 

(5) 

and the parameters A and e are chosen variationally through a stationary phase approxima- 
tion, then the hard constraint u(s)^ = 1 is replaced by the thermally averaged constraint 
(u(s)^) = 1. This result had been conjectured by Lagowski et al. |jl3[- We also established 
that imposing the constraint (u(s)^) = 1 for every point on the chain requires two param- 
eters; specifically, the parameter e is required to suppress the fluctuations at the ends of 
the chain |[TT| , [T3[| . Furthermore, if one uses only the variational parameter A, as suggested 



elsewhere [jT4[, then one obtains a Gaussian expression for G{R; L) which is wrong. In more 
recent work we showed that the same strategy, i.e., replacing the exact thermal weight 
by Eq. (|^) and evaluating A and e variationally, produces excellent quantitative results for 
the elastic response of a semifiexible chain under tension. Notice that the optimal values of 
A and e depend on the property of interest. 

Following these earlier works ||TT| , p!5[| we shall calculate G(-R; L) by replacing the true 
thermal weight in Eq. by \['Afi7'[u(s)] in Eq. and use a stationary phase approximation 
to set A and e. The equation for G{R] L) with the weight given by \E'mf[u(s)] is 



G{R] L) = n [dX [de [ V[u{s)] exp J -k ■ R - u{s)ds 

J-ioo (ZTlZr J J J Jo 
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where Vt represents a normalization constant. The functional integral over u(s) in Eq. (^) is 
done by replacing u by v = u — k/2A. The resulting path integral corresponds to a harmonic 
oscillator that makes a transition from u(0) — k/2A to u(L) — k/2A in imaginary " time" 
L. Using the standard result for the harmonic oscillator propagator the distribution 
function G{R] L) becomes 



G{R] L)= J dxj deexp(-F(A,e)) 



(7) 
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where 



/'AL-2e\ , f smh(ujL)\ , , ^, 



up 



with 



„ = _ + ^coth(-) (9) 

e La; , f^^LX 1 , . 

^_+^tanh(-)+^^ (10) 

e La; , /a;L\ . , 

,__ + ^,anh(-) (11) 



and 



For the case of large L the function F(A, e) may be written in the somewhat more manageable 
form 

F(A, e) = L/(A) + ^ log(2e + Z^^) - 2e - ^ log(/p^) + O + const. (13) 

where /(A) = (3/2)a; — A(l — r^) and r = R/L. Thus, as L becomes larger the function 
F(A, e) becomes more sharply peaked and the stationary phase approximation becomes more 
accurate. For leading order in L the stationarity condition for A is /'(A) = 0, which gives 



4 VI - 



(14) 



and we find that e does not show any significant, leading order variation with r. Substi- 
tuting the stationary values of A and e into Eq. (|13]) for F(A, e) yields the desired simple, 
approximate expression for G{R; L), 

N f 9L \ 

where iV is a normalization constant and we have again taken the leading orders in L in the 
exponential and in the prefactor. In our earlier work [11| we showed that the stationary phase 
approximation reduces the persistence length from Ip, which obtains in an exact treatment of 
the thermal weight in Eq. (^, to Imf = (2/3)/p. Identifying the effective persistence length 
hiF with the measured persistence length of the polymer produces a simple approximate 
expression for the radial probability density of a semifiexible chain in three dimensions, 

/ X 4:TTNr^ / 3t 1 \ 
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where t = L/Imf- The normahzation constant is determined by the requirement that 







1 

P{r-t)dr = l. (17) 



The integral can be evaluated by the substitution r = x/yl + x^ to yield 

4a^/^e" 

~ 773/2(4 + 12a-i + 15a-2) V ) 

where a = 3t/4. 

The distribution function, P(r; t), vanishes as as r ^ and also vanishes at r = 1. The 
peak of the distribution function occurs at r^ax = [{v+Vv^ + 14)/7]^/^, where r] = (5/2) — a. 
As expected, when t — 0, then Vmax — > 1 and P(r, t) — > 5(r — 1). In Fig. (1) we plot P(r; t) 
for the five values of t for which WF [|lOl have presented simulation data. For comparison. 



the results of the analytic expressions obtained by WF [|T0| are also presented. For very stiff 



chains {t = 0.5 in Fig. (1)) the WF theory gives a more accurate estimate of the peak [17 . 
Nevertheless, over the range of t examined by WF |jlO[ our simple expression in Eq. (pIBf ) 
reproduces the data quite accurately. 

In order to further assess the validity of our theory we have also calculated the second 



and fourth moments using the radial distribution function in ([161) . The moments are given 
by 

/i„ = (r*^) = r^'Pir; t)dr. (19) 
Jo 

The same substitution that was used to evaluate the integral in Eq. (|T7|) can be used to find 
the moments, yielding 

/3\ 4a + 10 

= (2) 4„^ + 12„ + 15 (^°' 

15 

- Aa^ + 12« + 15 ^^^^ 

for the second and fourth moments. In the limit of a 00 (random coil limit) ^2 ~^ 3 /2a 
and /i4 — i> 15/4a^, both of which coincide with results using exact moments [T^. Similarly 
both /X2 and ^4 tend to unity in the limit of a ^ 1 (hard rod limit), which is once again the 
exact result. 

The theory outlined here provides surprisingly accurate results for P{r;t). The correct 
limiting behavior is obtained for P(r; t) and for the first two moments. Furthermore, our 
theory can be systematically applied to other problems involving semifiexible chains. Thus 
it appears that the meanfield-like approach used here may be useful for treating a wide 
variety of problems in which the rigorous enforcement of the constraint u(s)^ = 1 is difficult 
to enforce. 

We wish to thank Dr. G. H. Weiss for useful suggestions and a referee for pointing out 
a misprint in an earlier version of this paper. 
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Figure Caption 



Fig. 1: Comparison of P{r,t) obtained with Monte Carlo simulations (represented by sym- 
bols) for t = 10, 5, 2. 1, and 0.5 and with analytic theories (represented by curves) for the 
same values of t. The curves and symbols are arranged so that the largest value of t is the 
left most and the smallest value of t is the right most. The dark lines corresponds to Eq. 
(P!B|) and the light lines are based on the approximate theory of Wilhelm and Frey ||10[ |. 
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